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ABSTRACT 

This  paper  reports  on  analytical  work,  as  well  as  experimental  testing,  that 
were  accomplished  at  the  Ames  Research  Center  and  at  the  Marshall  Space  Flight 
Center  to  examine  acoustic  wave  propagating  and  the  ability  to  detect  intrinsic  faults 
in  sandwich  type  composite  structures.  Sandwich  type  composites  are  being  studied 
for  use  in  NASAs  new  heavy  lift  launch  vehicle  and  flaw  detection  is  crucial  for 
safety  and  for  failure  prognostics.  The  work  reported  on  in  this  paper  involved  both 
the  theoretical  modeling  as  well  as  comparison  with  empirical  testing  needed  to 
answer  the  question  of  feasibility  for  reliable,  and  accurate,  structural  health 
monitoring  (SHM)  in  the  composite  structure  of  interest.  The  analytical  model  of  the 
transient  wave  propagation  and  scattering  based  on  the  Mindlin  plate  theory  was 
developed.  A  scattered  transient  field  properties  are  calculated  theoretically  using 
this  model  and  numerically  using  the  finite  element  model  for  acoustic  waves 
generated  by  an  acoustic -patch  actuator.  It  is  shown  that  theoretical  results  are  in 
agreement  with  the  results  of  numerical  simulations  and  with  experimental  results. 

Introduction 

Composite  sandwich  panels  (CSP)  are  replacing  metals  in  engineering 
applications  due  to  their  high  strength  and  low  weight  properties,  as  well  as  their 
flexible  production  capabilities  [1].  Their  stiffness-to- weight  ratios  and  the  lack  of 
corrosion  lead  to  the  wide  usage  of  them  in  the  naval  ships,  airplanes,  bodies  for 
vehicles  and  trains,  etc.  They  are  very  prospective  for  space  applications  where 
light  flight  structures  immediately  lead  to  higher  payloads.  CSP,  to  which  we 
restrict  our  consideration  here,  are  complex  structural  materials  made  from  stiff 
facesheets  and  soft  cores  and  are  usually  referred  to  as  sandwiches  with  soft  cores. 
Sandwich  panels  have  one  basic  drawback  in  common:  they  all  share  failure  mecha¬ 
nisms  and  structural  failures  occur  internally,  out  of  the  view  of  normal  visual 
means  of  inspection.  Internal  de-bonds,  de-laminations,  cracks,  and/or  buckles  are 
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typical  local  failure  modes,  decrease  the  strength  of  composite  material  sandwich 
structures,  and  may  lead  to  serious  damages  and  even  to  structure  disintegration. 

Many  methods  are  used  for  investigating  acoustic  wave  propagation  and 
scattering  in  the  framework  of  structural  health  monitoring  (SHM).  We  consider 
that  the  waves  are  excited  by  piezoelectric  transducers  (PZT),  which  act  both  as 
actuators  and  sensors.  In  this  paper  we  use  plate  theory  for  analytical  investigation 
of  wave  propagation  in  sandwich  panels,  discuss  parameters,  governing  equations, 
and  the  properties  of  the  basic  sandwich  panel.  High-fidelity  simulations  are 
employed  to  compare  the  results  of  the  simulation  with  analytical  results.  Some 
comparison  with  experimental  study  is  made. 


Plate  theory  for  sandwich  structures. 


In  a  Mindlin  plate  theory  the  displacements  of  the  plate  in  the  transverse,  radial,  and 
tangential  direction  components  are  expressed  as  follows  [2-5] 

w  =  w(r,  0,  t),u  =  z  Yr  ( r >  ^  t),v  =  z  Ye  (r>  0, 


where  z  is  the  coordinate  defining  points  across  the  thickness  of  the  plate  ( z  =  0  is 
the  neutral  plane)  and  w  is  the  out-of-plane  displacement,  Yr  and  Ye  are  the 
rotations  of  vertical  lines  perpendicular  to  the  mid-plane.  The  governing  equations 
for  the  symmetric  honeycomb  panels  in  terms  of  moments  M  and  shear  forces  Q 
can  be  presented  as 
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where  p  =  V3_  I'*'  p.dz,  is  the  mass  density  per  unit  area  of  the  plate,  index  k 

corresponds  to  the  material  layer,  pk  is  the  density,  /  =  V3_  \bk  p,z2dz  is  the 

mass  moment  of  inertia.  Each  layer  in  the  sandwich  panel  is  bounded  by  the  coor¬ 
dinates  ah  and  bk  in  the  thickness  direction  as  shown  in  the  Fig.  1.  The  stress 
resultants  in  terms  of  moments  Mrr,  M  ee ,  and  Mre,  along  with  shear  forces  Qr 
and  Qe  can  be  related  to  the  transverse  displacements  and  rotations  as  follows: 
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where  D  -  Eft3f /6  + Ect3c /\2  + +tc)2  jA-  is  the  flexural  stiffness,  vis  the 
Poisson  ratio,  which  for  the  sake  of  simplicity  is  taken  as  equal  for  each  layer,  Ef,  Ec 


-  Young's  modulus  of  the  facesheet  and  the  core,  correspondingly,  tf,  tc  -  thicknesses 
of  the  facesheet  and  the  core  layers,  G  is  the  shear  stiffness  of  the  plate. 

Wave  propagation.  The  general  solution  of  the  acoustic  waves  propagation 
with  cycling  frequency  co  is  (w,y/)  =  /te[(W,vF)exp(-haO]  [2-7],  where  and 

throughout  this  paper,  Re{.)  denotes  the  real  part  of  the  quantity  in  parentheses, 
y/=(  yr.  Ye)-  For  isotropic  sandwhich  layers  wavenumbers  have  three  brunches: 


(5) 


where 


The  flexural  wave  corresponds  to  the  real  (0,k  in  whole  co  domain.  The  second  (and 
third)  dilatation  branch  of  (k,  co)  dependence  become  real  starting  from  the  cutoff 
frequency. 

Model  of  the  scatterer.  The  scattering  problem  for  Mindlin  plate  theory  was 
considered  in  several  articles.  Pao  and  Chao  [6]  considered  scattering  of  a  flexural 
wave  at  a  circular  hole.  The  circular  inhomogeneity  was  considered  in  [4]  and  the 
graduate  damaged  model  -  in  [7], 

Let  us  consider  the  following  scattering  problem  (see  Fig.  lb).  The  circular 
delamination  (debond)  is  located  at  the  center  of  the  sandwich  composite  panel 
(x=y=0).  The  source  of  the  cylindrical  wave  generated  by  PZT  sensor  is  located  at  a 
distance  b  from  the  center  of  the  circular  delamination  (x=  -  b,  y=0).  In  this  paper 
we  extend  the  damage  model  of  the  defect  in  the  CSPs  in  such  a  way  that  it  allows 
us  to  differentiate  between  the  delamination  and  debonding  faults.  The  sketch  of  the 
damage  model  is  shown  in  Fig.  la.  It  consists  of  two  vertically  stack  circular 
cylinders  with  integral  properties  depending  on  the  location  of  the  fault.  If  the  fault 
is  located  in  the  facesheet  the  model  represents  a  delamination.  If  the  fault  is 
located  at  the  interface  between  the  facesheet  and  the  soft  core  the  model  represents 
a  debond.  To  find  the  parameters  of  the  fault  we  will  use  formulas  for  a  non— 
symmetrical  sandwich  panel  [1],  which  for  brevity  are  not  presented  here.  In  this 
case,  the  symmetrical  sandwich  part  of  the  panel  is  denoted  by  index  1.  The 
delaminated  region  consists  of  two  subregions  (see  Fig.  lb):  (i)  a  subregion  with  the 
sandwich  core  denoted  by  index  2  and  (ii)  the  delaminated  facesheet  subregion 
denoted  by  index  3  with  plate  parameters  D2,  h,  G2,  p2  and  Z)3,  /3,  G3,  p3, 
correspondingly.  If  we  introduce  wavenumbers  ky  for  the  delaminated  regions  we 
can  use  the  dispersion  relation  for  the  roots  (5)  with  parameters  corresponding  to 
indices  2  and  3  (index  1  belongs  to  the  symmetric  sandwich  panel  without  any 
damage).  As  a  result,  we  get  a  square  matrix  of  k(j  where  the  first  index  is  a 
number  of  a  root  (5)  and  the  second  one  is  a  number  of  a  media  (/=  1,2,3  ).  For  the 
self-consistent  behavior  of  the  delaminated  region  we  impose  boundary  conditions 
between  the  delaminated  region  and  the  healthy  annular  sandwich  beam.  These 
conditions  include  the  continuity  of  the  plate  displacements,  the  equality  of  the 
slopes,  moments  and  shear  forces  at  r=a: 


(6) 


r=a 


b2la: 


x= — by  j  ! 


scatterer 


Figure  1.  Sketch  of  the  sandwich  structure  with  delamination  fault.  Domain  1  (healthy  panel)  is 
connected  with  two  domains  2  and  3  through  the  boundary  of  the  fault  -  a,  Geometry  of  the 
scattering  of  the  incident  wave  -  b  . 


where  moments  and  shear  forces  are  described  by  expressions  (4),  with 
corresponding  indices. 

Incident  wave.  We  will  consider  a  circular-patch  actuator  on  the  Mindlin 
plate  generated  by  a  surface  traction  plate  waves  in  the  form  derived  in  [2,5],  For 
the  sake  of  brevity  let  us  consider  only  w'n  component  which  for  the  actuation 
circular  force  leads  to  the  incident  wave  for  r<b 

win(r,  03)  =  ^e„(-l)"H„(^)/„a„r)cos(«0),  (7) 

n= 0,  i=l 

where  s,  are  a  certain  constants  depending  on  wavenumbers  k,  and  Hankel  functions 
H0(k\r  +  b\)  are  expanded  into  a  Fourier  series  and  £,,  is  the  Neumann  factor,  £,=  1 , 
£,i=2  (n=  1,2,3,...)  [8], 

The  solution  of  scattering  waves  from  the  fault.  If  we  consider  that 
equations  (l)-(4)  are  true  for  each  domain:  outside  the  damage  and  inside  of  it 
(/'=  1,2,3),  the  solution  is  separated  into  three  domains,  where  the  delaminated 
region,  in  turn,  consists  of  two  domains.  As  a  result,  the  general  solution  can  be 
represented  for  r  >  a  by  incident  and  scattered  waves.  We  seek  the  solution  in  the 
form  of  a  linear  combination  of  the  scattered  and  incident  waves.  The  functions  we 
have  used  for  the  solution  should  be  finite  at  r  -  0  and  vanish  at  r  — >  °°.  As  a 
result,  the  general  solution  for  the  domains  i  =  1,2,3  is 

wl  =  w‘nn  +  alnHn(kur)  +  a2nHn(k21r)  \r>a 

W2  =  «4  Jn  (V)  +  «5  nJ  n  ( Kir )  K<a  (8) 

y  =  a  4 n  Jn  (^13r)  +  Cl ' 5n  Jn  (k23r)  l<a 

Here  wavenumbers  %  are  calculated  for  three  media  and  the  first  index  i 
corresponds  to  a  root,  and  the  second  index  j  -  to  a  medium  (expressions  for  y/r,  ye 
are  not  presentet  here  for  making  writing  concise).  The  boundary  conditions  (6), 
lead  to  a  system  of  linear  algebraic  equations  for  each  value  of  n 

\Un=bn  (9) 

The  matrix  An  and  vector  bn  and  complete  list  of  the  matrix  elements  can  be 
calculated  by  symbolic  calculus  and  is  not  presented  here. 

As  a  result  we  obtain  a  system  of  linear  algebraic  equations  relative  to 


w(r,t,6 )  =  Re'y'g  l&*cos (n0)< 

n= 0 


constants  un  =  (a]n ,  a2n ,  a3n ,  aAn ,  aAn ,  a5n ,  a5n ,  a6n ,  a6n  )T ,  which  can  be  solved  nume¬ 
rically.  We  find  that  the  first  several  main  modes  determine  the  spectrum  of  the 
scattered  wave,  and  truncation  of  n  at  n=30-50  in  eqautions  (9)  gives  practically  the 
same  result  at  any  desired  wave  frequencies.  According  to  the  expression  of  a\n. 


«!„(©)  = 


det  Aln  (ky  ((0)a) 
det  An  (kj(a)a) 


a2n  (to)  = 
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the  elastic  field  around  the  circular  debond  is  computed  and  the  coefficients 
aln,a2n  determine  displacements  outside  the  debond  region.  In  expressions  (10) 
matrix  Au(kija)  is  obtained  by  substituting  the  first  column  in  matrix  An  (ktja) 
by  the  left  hand  side  column,  and  matrix  A2n(kija)  -  by  substituting  the  second 

column  by  the  right  hand  side  column.  As  a  result,  for  out  of  plane  solution  for 
r  >  a  ,  the  formula  can  be  presented  as 


w(r, 6, co)  =  Re 2>os(n0)[w; (r)  +  ^ (kn{co)r)  +  a2n (< <o)Hn (k2l  {co)r))e~iM  (11) 

n= 0 

The  first  term  in  (11)  corresponds  to  the  incidend  wave,  second  -  to  the  flexural 
mode  with  a  wavenumber  kw  and  the  thrid  one  describes  the  flexural  wave  with  a 
wavenumber  &21  propagating  above  the  cutoff  frequency.  Below  the  cutoff 
frequency  (6)<  G\)  the  third  term  corresponds  to  the  evenecent  wave.  The  solution 
for  rotations  of  vertical  lines  perpendicular  to  the  mid-plane  if/,  and  y/e  can  be 
determined  in  a  similar  way. 

Transient  Solution.  To  study  the  transient  wave  propagation  we  consider 
that  the  plate  is  excited  by  a  pulse  of  the  load  stimulated  by  a  PZT  sensor.  The 
expression  for  the  wave  pulses  in  the  plane  (x,y)  may  be  derived  from  the 
steady- state  solution  in  the  frequency  domain  by  applying  the  Fourier  transform 
technique.  For  the  Fourier  spectrum  G(co)  of  the  signal  g(t)  the  final  solution  for  out 
of  plane  field  outside  the  debond  in  the  time  domain  will  be 


w(r,0,t)  =  Re — 
271 


|  G(co)w(co,  1 


We  consider  the  Hanning  type  actuation  signals  which  are  usually  used  for 
fault  detection  in  SHM.  Out  of  plane  displacemnet  (8)  is  an  analytical  solution  of 


Figure  2.  Wave  propagation  and  scattering  by  the  debond  region  located  at  the  center  of  the 
plate  (dashed  line  -  pristine  plate,  solid  black  line  -  plate  with  delamination:  v=-13.5cm  -a), 
a --5 cm  -  b),  a  -  1 3cm  -  d). 


scatteing  domain  -  a);  Scattering  from  the  debond  region  located  at  the  center.  Path  is  taken  along 
the  line  passing  through  the  centers  of  the  source  and  the  scatterer  -  b). 

the  scattering  problem  and  makes  it  possible  to  investigate  the  transient  pulse 
propagation,  and  to  compare  the  obtained  results  with  the  finite  element  simulation. 

In  Fig.  2  the  transient  pulse  propagation  determined  by  formulas  (8), (12)  is 
presented  for  different  points  of  the  line  passing  through  the  actuator  and  the  fault. 
The  dynamics  of  the  propagating  signal  for  pristine  plate  for  these  points  of  the 
structure  is  plotted  by  the  blue  dashed  line  and  the  w  displacement  of  the  debond 
plate  is  depicted  by  the  black  line.  We  can  see  that  for  the  pristine  plate  without  a 
fault,  the  wave  propagates  without  any  reflection  and  due  to  the  dispersion  we  have 
an  increasing  number  of  oscillations  (Fig.  3  a).  In  the  plate  with  a  fault  we  see  that 
after  scattering  the  propagating  signal  partially  propagates  back  and  continues  to 
propagate  forward  in  the  initial  direction  (Fig.  3  b).  The  transient  solution  of  the 
propagating  wave  as  a  result  of  calculation  of  the  expression  (7)  shows  that  the 
form  of  the  signal  is  different  for  each  specific  point  of  the  plate.  While  analyzing 
the  structure  of  the  transient  signal  we  can  see  two  main  features.  The  signal  for  the 
points  located  sufficiently  far  from  the  fault  consists  of  the  direct  signal  and  the 
reflected  one  of  the  similar  form  propagating  back.  Signals  passing  along  the 
straight  line  connecting  the  actuator  with  the  fault  have  counter  motion  as  it  is  seen 
from  Figures  2  and  3.  The  measurement  of  the  back-propagating  signal,  which  can 
be  detected  by  the  pulsed  echo  technique,  and  the  forward-propagating  one  detected 
by  the  pitch-catch  method,  can  retrieve  information  about  the  fault  since  they  carry 
this  information  while  propagating  through  the  plate. 

When  transient  waves  are  used  to  detect  delamination  in  the  sandwich 
plate,  the  total  solution  can  be  written  as  a  superposition  of  the  incident  and 
scattering  waves  in  the  time  domain  containing  n  + 1  different  modes  (7). 


Figure  4.  FE  modeling  of  wave  propagation  from  the  center  and  scattering  from  the  debond  region 
(view  in  x,y  plane,  damage  is  located  at  x=14cm  fro  the  center  (0,0))  at  £=0.45xl0~4s  -  a), 
f=0.9x!0'4i  -b,  f=1.35xl0'4s  -  c)). 


We  compare  the  analytical  results  with  the  corresponding  results  obtained 
by  the  Finite  Element  simulation.  The  FE  modeling  for  2D  Mindlin  plate  is 
presented  in  Fig.  4.  and  it  fits  well  with  the  theoretical  approach.  The  main 
difference  between  these  signals  is  that  the  theoretical  results  are  valid  for  the 
infinite  plate,  and  the  FE  2D  Mindlin  plate  model  gives  the  result  which  takes  into 
account  reflections  from  the  boundaries.  The  comparison  of  the  3D  modeling  of 
sandwich  composite  panels  with  the  theoretical  result  is  considered  in  the  last 
section.  The  propagation  of  the  wave  originated  at  the  center  of  the  plate  (0,0)  and 
scattered  from  the  debond  region  located  at  the  distance  b=0.14m  apart  is  presented 
in  Figure  4.  You  can  see  that  the  cylindrical  out-of-  plane  wave  propagates  without 
any  reflections  (Figure  4  a),  and  the  scattered  from  the  debond  region  wave  excites 
the  vibration  of  the  debond  generating  a  pulsed  echo  signal  and  a  changing  pitch 
catch  signal. 

Experimental,  Analytical  and  FE  comparisons.  An  experiment  in  Lamb 
wave  propagation  in  a  honeycomb  sandwich  panel  was  done  by  Metis  Design  Inc. 
in  collaboration  with  ARC  NASA.  The  sandwich  panel  fabricated  for  this  test 
consisted  of  two  84-mil  thick  cross-ply  carbon  fiber  composite  laminates  (bonded 
to  a  1”-  thick  aluminum  honeycomb  core).  The  size  of  the  panel  was  lftxlft.  In  the 
experiment,  PZT  sensors  located  on  the  facesheet  of  the  honeycomb  panel  were 
used  to  determine  the  deformation  at  a  different  location  (See  disks  located  on  the 
panel  which  correspond  to  real  PZT  sensors  in  the  experiment).  The  results  of  the 
comparison  voltage  on  these  PZT  sensors  are  presented  in  Figure  5  b,  c. 

The  obtained  experimental  time  domain  results  are  compared  with  the 
corresponding  results  of  the  FE  simulations  in  Abaqus.  The  FE  analysis  makes  it 
possible  to  find  the  amplitude  of  all  displacements  (, u,v,w )  in  3D  simulations  as  well 
as  the  stresses  and  strains  under  an  external  excitation  depending  on  the  frequency 
of  the  excitation.  The  FE  model  was  built  as  an  exact  copy  of  the  experimental 
panel.  It  was  meshed  to  define  the  solution's  resolution  and  satisfy  the  tolerance 
conditions,  which  make  it  possible  to  investigate  the  acoustic  wave  propagation 
through  a  layered  sandwich  plate.  The  discretization  was  chosen  to  determine  the 
wavelength  with  a  high  resolution  at  several  elements  per  wavelength  limits.  Finite 
element  3D  modeling  of  real  honeycomb  structure  with  laminate  facesheets  is 
performed  in  order  to  validate  the  approximating  Mindlin  plate  theory  analytical 
approaches  and  to  get  more  realistic  results  in  application  of  the  homogenized 
structure  model  for  SHM.  The  finite  element  simulation  of  Abaqus  Explicit  has 


Time,  sec  x  10’5  Time,  sec  xio'5 

Figur  5.  View  of  the  panel  FE  model  used  for  experimental  investigation  of  wave 
propagation  and  scattering  -  a);  Sensor  signal  in  mV  on  the  sensors  located  at  a  distance  of  25cm 
apart  from  each  other  for  a  sandwich  composite  panel:  x=-12.5  cm  -b  and  12.5cm  -  c.  The  red 
line  -  experimental  data,  black  line  -  FE  simulations  and  the  green  line  -  fitted  Mindlin  plate  theory 
results. 


been  used  to  determine  the  PZT  actuation  of  the  sandwich  plate.  The  small  debond 
fault  was  introduced  at  the  center  of  the  plate  (Fig.  5  a).  For  the  sandwich  structure 
we  performed  the  simulation  of  the  several  possible  layups  we  come  across  in 
developing  SHM  methods.  The  actuator  and  the  sensor  were  simulated  as  PZT  discs 
(14mm  and  7mm  in  diameter  and  2mm  in  thickness).  The  excitation  source  was 
positioned  at  a  distance  b  from  the  damaged  area  and  the  receiving  sensors  are 
located  in  (x,  y)  plane  as  it  is  seen  from  the  plots. 

In  the  Mindlin  plate  theory  we  do  not  consider  coupling  with  piezo  sensors 
and  that  is  why  we  cannot  plot  a  real  electric  signal  on  this  plot.  We  used 
comparing  relative  amplitude  and  phase  velocity  to  show  how  the  Mindlin  plate 
theory  results  relate  to  a  real  structure.  The  comparison  of  out  of  plane  displacement 
dynamics  shows  that  the  obtained  results  are  a  sufficiently  good  fit  to  analytical 
results.  The  simulations  are  done  for  modeled  structures  with  lin  core  and  84mil 
face  sheets  in  thickness. 

It  should  be  mentioned  that  the  results  of  modeling  wave  propagation  and 
scattering  in  sandwich  composite  plates  fit  well  to  the  theoretical  ones  only  for 
small  time  intervals.  For  a  longer  time  interval,  the  experimental  result  differs 
significantly  from  the  theoretical  approach.  However,  if  in  the  case  of  FE  modeling 
we  can  fit  experimental  and  numerical  signals  for  a  longer  time,  here  the  theoretical 
approach  is  proposed  for  the  infinite  plate  only. 

Conclusion 

This  paper  addresses  the  simulation  of  the  guided  wave  propagation  and  scattering 
in  sandwich  structures  with  the  emphasis  given  to  the  properties  which  can  be  used 
for  SHM.  The  analytical  investigation  of  the  plate  wave  using  the  Mindlin  plate 
theory  and  the  numerical  simulations  shows  the  main  features  we  come  across  when 
developing  real  SHM  methods.  An  analytical  study  is  carried  out  to  find  the 
solution  for  transient  wave  propagation  and  scattering  from  the  delaminated  and  the 
debond  area.  The  obtained  analytical  solutions  are  compared  to  the  FE  analysis  as 
well  as  the  experimental  data. 
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